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Abstract. The Chandra X-ray Observatory is finding a large number of cavities in the X-ray emitting intra-cluster medium 
which often coincide with the lobes of the central radio galaxy. We propose high-resolution Sunyaev-Zel'dovich (SZ) observa- 
tions in order to infer the still unknown dynamically dominating component of the radio plasma bubbles. This work calculates 
the thermal and relativistic SZ emission of different compositions of these plasma bubbles while simultaneously allowing for 
the cluster's kinetic SZ effect. As examples, we present simulations of an Atacama Large Millimeter Array (ALMA) observa- 
tion and of a Green Bank Telescope ( GBT) observation of the cores of the Perseus cluster and Abell 2052. We predict a 5cr 
detection of the southern radio bubble of Perseus in a few hours with the GBT and ALMA while assuming a relativistic electron 
population within the bubble. In Abell 2052, a similar detection would require a few ten hours with either telescope, the longer 
exposures mainly being the result of the higher redshift and the lower central temperature of this cluster. Future high-sensitivity 
multi-frequency SZ observations will be able to infer the energy spectrum of the dynamically dominating electron population 
in order to measure its temperature or spectral characteristics. This knowledge can yield indirect indications for an underlying 
radio jet model. 
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1. Introduction 

The Chandra X-ray Observatory is detecting numerous X-ray 
cavities in clusters of galaxies confirming pioneering detec- 
tions of the ROSAT satellite. Prominent examples include the 
Perseus cluster (Bohringer et al. 1993; Fabian et al. 2000), 
the Cygnus-A cluster (Carilli et al. 1994), the Hydra-A clus- 
ter (McNamara et al. 2000), Abell 2597 (McNamara et al. 
2001), Abell 4059 (Huang & Sarazin 1998; Heinz et al. 2002), 
Abell 2199 (Fabian 2001), Abell 2052 (Blanton et al. 2001), 
the vicinity of M84 in the Virgo cluster (Finoguenov & Jones 
2001), the RBS797 cluster (Schindler et al. 2001), and the 
MKW3s cluster (Mazzotta et al. 2002). They are produced by 
the release of radio plasma from active galactic nuclei (AGN) 
which are typically hosted by a cD galaxy located at the cluster 
center and mainly reside within the cool cores of galaxy clus- 
ters. 

While radio synchrotron emission provides evidence for the 
existence of cosmic ray electrons (CRes) and magnetic fields, 
the detailed composition of the plasma bubble governing its dy- 
namics is still unknown. Minimum energy or equipartition esti- 
mates of the nonthermal pressure in the radio bubbles give val- 
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ues which are typically a factor of ten smaller than the pressures 
required to inflate and maintain the bubbles as determined from 
the surrounding X-ray gas (e.g., Blanton et al. 2001). This indi- 
cates that the standard minimum energy or equipartition radio 
arguments are missing the main component of the pressure and 
energy content of the radio lobes. Possibilities include mag- 
netic fields, cosmic ray proton (CRp) or CRe power-law dis- 
tributions, or very hot thermal gas. Solving this enigma would 
yield further insight into physical processes within cool cores 
(De Young 2003) as well as provide hints about the composi- 
tion of relativistic outflows of radio galaxies because plasma 
bubbles represent the relic fluid of jets (e.g. Celotti et al. 1998; 
Hirotani et al. 1998; Sikora & Madejski 2000). 

Additionally, some of the clusters exhibit cavities in the 
X-ray emitting intra-cluster medium (ICM) without detectable 
high frequency radio emission, for instance in Perseus, Abell 
2597, Abell 4059, and the MKW3s cluster. This category of 
X-ray cavities is also believed to be filled with radio plasma, 
but during the buoyant rise of the light radio plasma bubble in 
the cluster's potential (Gull & Northover 1973; Churazov et al. 
2000, 2001; Briiggen & Kaiser 2001) the resulting adiabatic 
expansion and synchrotron/inverse Compton losses dwindles 
the observable radio emitting electron population producing a 
so-called ghost cavity or radio ghost (EnBlin 1999). Possible 
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entrainment of the ICM into the plasma bubble and subsequent 
Coulomb heating by CRes generates further uncertainty of the 
composition of the ghost cavity. 

If the radio bubbles contain a significant amount of very 
hot thermal gas, this might be detected by X-ray observations. 
However, this is quite difficult (e.g., Blanton et al. 2003) due to 
the projected foreground and background cluster emission, and 
the fact that the X-ray emissivity is proportional to the square 
of the density. If most of the pressure in the radio bubbles were 
due to the very hot, low density thermal gas, it would have a 
very low X-ray emissivity. Observationally, there has been a 
claim by Mazzotta et al. (2002) to have seen hot X-ray emitting 
gas within the ghost cavity of the MKW3s cluster. Obviously, 
it might be more useful to observe the radio bubbles with a 
technique which was sensitive to thermal gas pressure, rather 
than density squared, as pressure is the quantity which is miss- 
ing. For this reason, in this paper we propose high resolution 
Sunyaev-Zel'dovich (SZ) radio observations of the radio bub- 
bles and radio ghosts in clusters, as the thermal SZ effect di- 
rectly measures the thermal electron pressure in the gas. The 
thermal SZ effect arises because photons of the cosmic mi- 
crowave background (CMB) experience inverse Compton col- 
lisions with thermal electrons of the hot plasma inside clus- 
ters of galaxies and are spectrally redistributed (e.g. Sunyaev & 
Zel'dovich 1972; Sunyaev & Zeldovich 1980; Rephaeli 1995a). 
The proposed measurement is able to infer directly the compo- 
sition of radio plasma bubbles and radio ghosts while indirectly 
obtaining indications for a specific underlying jet model. 

The article is organized as follows: After basic definitions 
concerning the thermal, kinetic, and relativistic SZ effect in 
Sect. 2, we introduce a toy model in Sect. 3 describing pro- 
jected maps of the SZ flux decrement with spherically symmet- 
ric radio plasma bubbles. The models for the cool core regions 
of the Perseus cluster and Abell 2052 are described in Sect. 4. 
Simulating an Atacama Large Millimeter Array (ALMA ) and 
a Green Bank Telescope (GBT) observation of both clusters 
in Sect. 5, we examine whether the plasma bubbles are de- 
tectable by the SZ flux decrement. Five physically different 
scenarios for the plasma composition of the bubbles are inves- 
tigated exemplarily using three characteristic SZ frequencies 
in Sects. 6 and 7. Finally, we discuss observing strategies for 
ALMA and GBT. Throughout the paper, we assume a ACDM 
cosmology and the Hubble parameter at the present time of 
H = 70/z 7 o km s" 1 Mpc" 1 . 

2. Sunyaev-Zel'dovich effect 

The SZ effect arises because CMB photons experience inverse 
Compton (IC) scattering off electrons of the dilute intra-cluster 
plasma (for a comprehensive review, see Birkinshaw 1999). At 
the angular position of galaxy clusters, the CMB spectrum is 
modulated as photons are redistributed from the low-frequency 
part of the spectrum below a characteristic crossover frequency 
v c to higher frequencies. For non-relativistic electron popula- 
tions, v c ^ 217 GHz, while this characteristic frequency shifts 
towards higher values for more relativistic scattering electrons. 

The relative change Si(x) in flux density as a function of 
dimensionless frequency x = hv/ikTcya) for a line-of-sight 



through a galaxy cluster is given by 
8i(x) = g(x)y g . ds [1 + 6(x, T e )] -h(x)w m& 

+ [j(x)-i( X )]T rel , (1) 

with the Planckian distribution function of the CMB 

x 3 

I(x) = i i(x) = ;' — — -, (2) 
e* - 1 

and i'o = 2(£7cmb) 3 I (he) 2 where 7cmb = 2.725 K, k, h and c 
denote the average CMB temperature, Boltzmann's constant, 
Planck's constant, and the speed of light, respectively. 

The first term in Eqn. (1) arises because of the thermal mo- 
tion of non-relativistic electrons (thermal SZ effect) and gives 
the spectral distortion 

(e*-l) 2 \ e- 4 -l / 

The amplitude of the thermal SZ effect is known as the thermal 
Comptonization parameter y gas that is defined as the line-of- 
sight integration of the temperature weighted thermal electron 
density from the observer to the last scattering surface of the 
CMB atredshiftz = 1100: 



= I dln e 

m e c l J 



;kT, 



(4) 



Here, cr T denotes the Thompson cross section, m e the electron 
rest mass, T e and n ejgas are electron temperature and thermal 
electron number density, respectively. For non-relativistic elec- 
trons the relativistic correction term is zero, 6(x, T e ) = 0, but for 
hot clusters even the thermal electrons have relativistic correc- 
tions, which will modify the thermal SZ effect (Wright 1979). 
These corrections have been calculated in the literature (see e.g. 
Rephaeli 1995b; Enfilin & Kaiser 2000; Dolgov et al. 2001; 
Itoh & Nozawa 2004), and can be used to measure the clus- 
ter temperature purely from SZ observations (e.g. Hansen et al. 
2002). 

The second term in Eqn. (1) describes an additional spectral 
distortion of the CMB spectrum due to the Doppler effect of 
the bulk motion of the cluster itself relative to the rest frame of 
the CMB. If the component of the cluster's peculiar velocity is 
projected along the line-of-sight, then the Doppler effect leads 
to a spectral distortion referred to as the kinetic SZ effect with 
the spectral signature 



h(x) = 



x 4 e x 



(5) 
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The amplitude of the kinetic SZ effect depends on the kinetic 
Comptonization parameter Wgas that is defined as 



Wa 



j dln e 



(6) 



where r gas is the Thomson optical depth, v gus is the average 
line-of-sight streaming velocity of the thermal gas, y6 gas = 
Vgas/c, and y6 gas < if the gas is approaching the observer. 

Finally, the third term in Eqn. (1) takes account of Compton 
scattering with relativistic electrons that exhibit an optical 
depth of 



Trel 



= £T T J"d/n e ,rei. 
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The flux scattered to other frequencies is i(x)T K \ while j{x)T m \ 
is the flux scattered from other frequencies to x = hvl(kTcMe)- 
It is worth noting, that in the limit of ultra-relativistic elec- 
trons and for x < 10, one can neglect the scattered flux, be- 
cause j(x) <sc i(x). In the following, we drop this approxima- 
tion and consider the general case. The scattered flux can be 
expressed in terms of the photon redistribution function for 
a mono-energetic electron distribution P(t; p), where the fre- 
quency of a scattered photon is shifted by a factor t: 

r*oo r*oc 

j(x)= dt dpMp)P(t;p)i(x/t). (8) 
Jo Jo 

For a given electron spectrum f e (p)dp with the normalized 
electron momentum p = /3 e y e and J dp f e (p) = 1, this redistri- 
bution function can be derived following the kinematic consid- 
erations of Wright (1979) of the IC scattering in the Thomson 
regime, where y e hv «: m e c 2 is valid. We use the compact for- 
mula for the photon redistribution function which was derived 
by EnBlin & Kaiser (2000): 

P ^p) = -^^[ 1 + ( 10 + 8 /> 2 + V)f + f 2 ] 
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The allowed range of frequency shifts is restricted to 
I ln(f)| < 2 arcsinh(/?) , 



(9) 



(10) 



and thus P(t; p) — for | ln(f)| > 2 arcsinh(p). Similar expres- 
sions for the photon redistribution function using different vari- 
ables can be found in the literature (Rephaeli 1995b; EnBlin & 
Biermann 1998; Sazonov & Sunyaev 2000). 

The spectral distortions owing to the relativistic SZ effect 
can be rewritten to include a relativistic Comptonization pa- 
rameter y, 



where 

y 

kf e 
g(x) 

P(kf e ) 



[j(x) - /(x)]T re i = g(x)y, 



= -^r- I dln e 



kf e , 



«e 

m,c 2 



i(x)]p(kT e ), 
m e c 2 J dl n e 
Wt) ~ ] dln e kT e 



(11) 

(12) 

(13) 
(14) 

(15) 



We introduced the normalized pseudo-thermal beta-parameter 
fi(kT e ) and the pseudo-temperature kT e which are both equal 
to its thermodynamic analog in the case of a thermal electron 
distribution. If the CRe population is described by the power- 
law distribution (29), the CRe pressure is given by 

PCRZ = -|- d P f(p)&p (16) 

j Jo 



n C R e m e c 2 (a - 1) 



la -2 3 - a\ 



(17) 



where f3 e = v/c = p/ y/l + p 2 is the dimensionless velocity 
of the electron, S q (a, b) denotes the incomplete Beta function 
(Abramowitz & Stegun 1965). In this case, the normalization 
of the CRe distribution function f(p)dp is determined by the 
CRe number density, «cRe = / dpf(p). Here, we introduced 
the abbreviation 



[F(p)C =F( Pl )-F( P2 ) 



(18) 



in order to account for the lower and upper cutoff p\ and p 2 of 
the CRe population. 

3. Toy model for plasma bubbles 

This section adopts an analytical formalism to describe buoyant 
plasma bubbles which was developed for the analysis of X-ray 
and radio emission by EnBlin & Heinz (2002). After a phase 
of supersonic propagation of the radio plasma into the ambient 
ICM, the radio lobes quickly reach pressure equilibrium with 
the surrounding medium once the AGN activity has terminated. 
During this stage, the bubble rises with constant velocity gov- 
erned by the balance of buoyancy and drag forces while the 
volume of the bubble expands adiabatically. Meanwhile, the 
surrounding gas is approximately in hydrostatic equilibrium 
with the underlying dark matter potential. Synchrotron, inverse 
Compton, and adiabatic losses diminish the observable radio 
emitting electron population within the plasma bubble produc- 
ing a so-called ghost cavity. 

As an analytically feasible toy model, we assume spherical 
geometry of the plasma bubble and adopt the general n-fold fi- 
profile for the electron pressure of the ICM which might find 
application for cool-core clusters: 



P e (r) = n e (r)kT e (r) = g P i 1 + |— j 



-3/W2 



(19) 



The origin of our coordinate system coincides with the cluster 
center while the x\- and X2-axes define the image plane, and the 
z-axis the line-of-sight to the observer. We choose the direction 
of the xi-axis such that the bubble center is located in the x\-z 
plane at r c = (r c cos 8, 0, r c sin 6). Its projected distance from 
the cluster center amounts to R c = pr c with p = cos 6, while 
the radius of the bubble is denoted by r^. For an unperturbed 
line-of-sight which is not intersecting the bubble, the observed 
thermal Comptonization parameter y c \(xi,X2) of the cluster is 
given by 



y c \{x u x 2 ) = J^: 



1 + 



r 2 . 

y,' 



-(3ft,,-l)/2 



■yb g , 



(20) 



where y,- = o- T (m e c 2 y l Pir y j!B ( 1 , |) is the central ther- 
mal Compton parameter of the respective individual y8-profile 
and ybg is the background contribution to the Comptonization 
which we set to zero in our analysis. In the case of a line-of- 
sight intersecting the surface of the bubble, the two intersection 
points are (xi , X2, z±) with 



P2 



Z± = r c ^l-p 2 ± -<\r\-x\- (x l - r c fi) 2 



(21) 
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The thermal Comptonization parameter y(x\ , X2) for the area 
covered by the bubble is given by 



>' b (xi , x 2 ) = y d (xi , x 2 ) - ^ yt 
sgn(z) 



1 + 



r 2 . 



-(3/? w -l)/2 



-I, 



(22) 



where I q = 1} q {a, b)/S(a, b) denotes the regularized Beta func- 
tion and q y j(z) = z 2 lir 2 yi + x\ + x\ + z 2 ). 

The amplitude of the kinetic SZ effect is proportional to 
the line-of-sight integrated electron density for which we also 
assume a general n-fold /3-profile: 



-3jS„,,/2 



(23) 



To avoid confusion, we adopt the notation r H .j and fi w j for the 
usual core radii and the /3-parameters in analogy to the ther- 
mal SZ effect. For an unperturbed line-of-sight which is not 
intersecting the bubble, the observed kinetic Comptonization 
parameter w c i(x\, x 2 ) of the cluster is given by 

-(3jS„,,-l)/2 



W c l(Xl,X 2 ) = ^\ 



1 + 



(24) 



where w,- = /3 g . d>i cr T njr w j!B ( 3/3 " 2 ' 1 , jj is the central kinetic 
Compton parameter of the respective individual /^-profile. The 
kinetic Comptonization parameter w(x\,x 2 ) for the area cov- 
ered by the bubble is obtained by analogy to the previous case: 

N ( v 2 _l v . 2 \-(3fe-l)/2 

w h (x u x 2 ) = w c \(xi,x 2 )-^ j w i 
sgn(z) 



i=i 



1 + 



2 

x\ + x 2 



2 

2 11JI 



-It 



|i 3Av,-l j 



2 2 

: , v-2 . ,2n 



(25) 



where g vv>i (z) = Z /(r^ + xf + 4 + z 2 )- 

4. Plasma bubbles of Perseus and Abell 2052 

Two of the most prominent examples of radio plasma bubbles 
in nearby galaxy clusters can be observed within the cool core 
regions of the Perseus cluster (redshift zp erS eus = 0.0179) and 
Abell 2052 (z A 2052 = 0.0348). Their proximity makes both 
clusters suitable targets for plasma bubble observations. Both 
clusters each host two bubbles which reflect the relic plasma 
of a past cycle of jet activity in the cD galaxy at the cluster 
center. At the current stage, the two radio lobes are in approx- 
imate pressure equilibrium with the surrounding medium and 
rise with a velocity governed by the balance of buoyancy and 
drag forces while the volume of the bubbles expands adiabati- 
cally. 

As in the previous section, the coordinate origin coincides 
with the cluster center and the new direction of the xi-axis 
points towards positive values of the relative right ascension. 
Assuming spherical symmetry of the plasma bubbles, their 
three dimensional position relative to the cluster center is de- 
generate because of projection effects. Together with the cluster 



center, the center of the bubbles form a plane which we assume 
to be perpendicular to the line-of-sight z. The azimuthal angle 
to the bubble center (p is measured from the xi-axis. 



4.1. Perseus 

Using deprojected electron density and temperature profiles de- 
rived from X-ray observations (Churazov et al. 2003), we ob- 
tain a pressure profile by fitting a double y0-profile according to 
Eqn. (19). In principle, we want the X-ray pressure profile in 
the absence of the radio bubbles. Since the observed X-ray sur- 
face brightness profile was derived assuming spherical symme- 
try including the region of the bubbles and doesn't extend into 
the very center of the cluster due to the AGN at the center and 
the XMM/Newton point spread function, our calculations some- 
what underestimate the SZ effect of the bubbles. Table 1 shows 
the individual parameters of the two plasma bubbles which are 
measured from the X-ray image of the central region of Perseus 
(Fabian et al. 2000). 

In our model, we adopt the peculiar velocity of the Perseus 
cluster of j} gas = -136 km s with respect to the rest frame of 
the CMB and approaching the observer (Hudson et al. 1997). 1 
The induced kinetic SZ effect gives rise to a small attenuation 
of the SZ decrement at our fiducial frequency. However, it also 
leads to an interesting effect at the crossover frequency v c 
217 GHz, and produces an enhancement of the SZ increment 
at higher frequencies (cf. Sect. 7). 

4.2. Abell 2052 

Despite the lower central pressure of Abell 2052 compared to 
Perseus, Abell 2052 lies at higher Galactic latitudes. Thus, SZ 
flux confusion with Galactic dust emission is negligible in this 
case, which might be an observational advantage. 

The electron density profile of Abell 2052 is obtained 
by deprojecting the X-ray surface brightness profile of Mohr 
et al. (1999) by means of the deprojection formula given in 
Appendix A of Pfrommer & Enfilin (2004). Since the X-ray 
surface brightness is represented by a double /?-model, the re- 
sulting density profile equals the square root of two single 
profiles added in quadrature. Refitting this electron density pro- 
file to match the profile defined by Eqn. 23 yields the parame- 
ters given in Table 1 . 

In order to model the temperature profiles T e (r) for 
Abell 2052, we applied the universal temperature profile for 



1 Using the fundamental plane, Hudson et al. (1997) measure 
the Malmquist bias corrected distance to early-type galaxies of the 
Perseus cluster. Comparing these distances to the mean of the indi- 
vidual galaxy redshifts, they infer the peculiar velocity of the galaxy 
cluster, L/ pcc = (—136 ± 307)kms~', where the uncertainty derives 
from the distance error added in quadrature with the cluster mean red- 
shift error. On the other hand, the predicted peculiar velocity derived 
from the IRAS redshift survey density field yields v pcc = +180 km s" 1 
where the velocity field was smoothed on a scale of 8 /T 1 Mpc which 
corresponds to the scale from which the clusters collapsed. 
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Table 1. Parameters of the deprojected electron pressure profiles which obey an n-fold /^-profile defined by Eqn. (19), the deprojected electron 
density profiles (Eqn. 23), and the individual parameters of the northern and southern plasma bubbles within the cool core regions of the Perseus 
cluster and Abell 2052, respectively. The parameters r c and r b correspond to the distance from the cluster center to the bubble center and the 
bubble radius, respectively. The azimuthal angle to the bubble center tp is measured from the axis defined by positive values of the relative right 
ascension while t] s denotes the angle of the normal vector of the mushroom-shaped southern bubble in A 2052 which we model as a half-sphere. 
In Perseus, 1" corresponds to 0.36 hz^ kpc and in A 2052, 1" corresponds to 0.69 kpc. 
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Cluster density profile: 
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cool core clusters proposed by Allen et al. (2001) to data taken 
from Blanton et al. (2001), 



W = T + (T] - T ) 



1 + 



F temp 



(26) 



where T = 1.31 keV, T x = 3.34 keV, r temp = 20.6 kpc, 
and 77 temp = 4.5. This equation matches the temperature profile 
well up to radii of ~ 0.3 r v ; r , which is sufficient for our pur- 
poses since we are especially interested in the core region of 
Abell 2052. Combining the electron density and temperature 
profiles yields the radial variation in the gas pressure, which 
we represent as a triple /^-profile (Table 1); for two of the com- 
ponents, the normalization is negative. As noted for Perseus, 
what we really need is the pressure profile in the absence of 
the radio bubbles. The observed X-ray surface brightness pro- 
file was derived assuming spherical symmetry including the re- 
gion of the bubbles. Probably as a result of this, the adopted 
pressure profile has more structure and a lower central pressure 
than would be true of the pressure profile in the absence of the 
radio bubbles. Thus, our calculations somewhat underestimate 
the SZ effect of the bubbles in Abell 2052. Since the southern 
bubble has a mushroom shape, we decided to model it using a 
half-sphere with the line-of-sight grazing the face. 

For our simulation, Abell 2052 is taken to be at rest in the 
CMB rest frame owing to the large uncertainties of the velocity 
determination with the fundamental plane method to measure 
the Malmquist bias corrected distance to early-type galaxies 
(Hudson et al. 2004; Hudson 2004). 2 



2 Hudson et al. (2004) have only 5 galaxies of Abell 2052 in 
their sample and infer a peculiar cluster velocity of w pcc = (-494 ± 
1062) km s _1 where the error is dominated by the error in the dis- 
tance but also includes redshift uncertainties and systematic effects 



5. Synthetic observations 

5.1. Atacama Large Millimeter Array 

In this section, only the general setup and the results for the 
synthetic ALMA observation of the cool core regions of Perseus 
and Abell 2052 are presented while the interested reader is re- 
ferred to Appendix A for further details. 

To simulate an ALMA observation of the central region 
of Perseus, we compute the frequency band-averaged SZ flux 
decrement <i5/) Vo around the frequency vo = 144 GHz which 
samples the extremum of the SZ flux decrement. The resulting 
simulated SZ flux decrement was convolved with a Gaussian of 
width <x Vo ^ 1 .2" to obtain the resolution of this ALMA config- 
uration. Since Perseus has more than 40 degrees of declination 
and is a Northern object, ALMA will be challenged to observe 
it, culminating at 25 degrees elevation. However, for a North- 
South elongated configuration of ALMA about the factor of 
(cos 65 °)-' ^ 2.4, the beam will be almost round. Assuming an 
ultra-relativistic electron population within the bubbles yields 
the SZ flux decrement as shown in Fig. 1 which is similar in 
morphology to the X-ray image. 

To investigate whether the plasma bubbles are detectable in 
the SZ flux decrement, we define the SZ flux contrast p: 



P = 



Ia-Ib 
Ia+Ib' 



(27) 



Here, I A and I B denote the mean SZ flux decrement of two 
equally sized solid angle elements within the field of view, one 
of which measures the SZ flux inside and the other one outside 
the bubble. Prior information about the angular position of the 



such as extinction. On the other hand, the predicted peculiar veloc- 
ity of Abell 2052 according the IRAS redshift survey density field is 
!V, = +208 km s -1 . 
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ALMA: Perseus cluster 
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Fig. 1. Synthetic ALMA observation of the cool core region of Perseus 
at the frequency band centered on 144 GHz. Shown is the simulated 
SZ flux decrement of radio plasma bubbles in units of mjy/arcmin 2 
assuming an ultra-relativistic electron population within the bubbles 
(scenario 1). The image is smoothed to the resolution of the ALMA 
compact core configuration (FWHM ^ 2.9"). The contour lines have 
a linear spacing of 0.16 mjy/arcmin". For comparison, the size of the 
ALMA field of view at this frequency is FWHM fov = 36" . 

X-ray cavities allows one to maximize the SZ flux contrast. 
Adopting Gaussian error propagation and introducing ALMA's 
sensitivity in terms of flux density, Oalma, yields the signal-to- 
noise for the 5cr detection of the bubble: 

s P Vaw^ (i 2 a - iQ 

H = ^r = , * 5. (28) 

9 2<talma Jl A + ^ B 

Here, M,eam is the number of statistically independent beams 
per flux averaged solid angle element (cf. Appendix A). We 
choose the field of view to be centered on the inner rim of the 
southern bubble such that equal solid angle elements fall in- 
side and outside the bubble, respectively. The simulated mean 
SZ flux decrements within the two solid angle elements of the 
primary beam require only an integration time of 5.1 hours in 
order to obtain a 5<x detection of the plasma bubble assum- 
ing an ultra-relativistic electron population within the bubble. 
Uncertainties in the amplitude of the kinetic SZ effect and the 
geometrical arrangement and shape of the bubbles may slightly 
modify this result. As a word of caution, this observation time 
allows a single bubble to be observed in a single pointing while 
it might be advisable to map the entire central region including 
both bubbles to get a clearer picture of the structure, and to be 
convinced that any holes seen in the SZ map at the radio bub- 



Table 2. Summary of predicted exposure times for the different com- 
binations of telescopes and clusters. / A and /b denote the mean SZ flux 
decrement of two equally sized solid angle elements within the field 
of view, one of which measures the SZ flux outside and the other one 
inside the bubble. 



Telescope: cluster 


Ia 


h 


exposure 




[mJy amin~ 2 ] 


[mJy amin~ 2 ] 


[hours] 


ALMA: Perseus 


13.25 


12.70 


5.1 


ALMA: Abell 2052 


3.930 


3.698 


38 


GBT: Perseus 


11.31 


10.85 


2.1 


GBT: Abell 2052 


3.272 


3.138 


31 



bles weren't just fluctuations also seen elsewhere in the cluster 
center away from the radio bubbles. 

The corresponding observation of Abell 2052 would re- 
quire an integration time of 38 hours in order to obtain a 5<x 
detection of the plasma bubble for the same plasma bubble con- 
tent, the longer exposures mainly being the result of the higher 
redshift and the lower central temperature of this cluster. 

5.2. Green Bank Telescope 

To simulate a GBT observation of the central region of both 
clusters, we compute the frequency band-averaged SZ flux 
decrement (6I) Vo in the frequency interval [86 GHz, 94 GHz]. 
The resulting simulated SZ flux decrement was convolved with 
a Gaussian of width <x 3.4" to obtain the resolution of the 
GBT 3 mm receiver. Assuming an ultra-relativistic electron 
population within the bubbles yields the SZ flux decrement as 
shown in Fig. 2. 

To investigate whether the plasma bubbles are detectable in 
the SZ flux decrement, we adopt the concept of SZ flux con- 
trast of the previous section. The sensitivity of the upcoming 
GBT Penn Array Receiver in terms of flux density is given by 
o~(jbt = 0.25 (Af) -1 ^ 2 mJy arcmin~ 2 , where At is the integration 
time in hours (Mason 2004). We choose the 32" x 32" field of 
view of the GBT to be centered on the inner rim of the south- 
ern bubble of Perseus as described above. Assuming an ultra- 
relativistic electron population within the bubble, we predict a 
5<x detection of the plasma bubble after an integration time of 
2. 1 hours owing to the better sensitivity of the bolometric re- 
ceivers on the GBT. This observation time assumes a proper 
foreground subtraction of the Galactic emission components. 

The corresponding observation of Abell 2052 would re- 
quire an integration time of 3 1 hours in order to obtain a 5cr 
detection of the plasma bubble for the same plasma bubble con- 
tent. Again, the integration times correspond to a single point- 
ing on the southern bubble while mapping the entire central 
region would take respectively longer. 

6. Composition study of plasma bubbles 

In the following, we study exemplarily five physically different 
scenarios of the composition of the plasma bubble which is as 
a whole in approximate pressure equilibrium with the ambient 
ICM. Although these scenarios might not be realized in nature 
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Fig. 2. Synthetic GBT observations of the cool core regions of Perseus (left panel) and Abell 2052 (right panel) at the frequency band centered 
on 90 GHz. Shown is the simulated SZ flux decrement of radio plasma bubbles in units of mjy/arcmin 2 assuming an ultra-relativistic electron 
population within the bubbles (scenario 1). The images are smoothed to the resolution the GBT 3 mm receiver (FWHM a 8.0")- The contour 
lines have a linear spacing of 0.16 mjy/arcmin 2 (left panel) and 0.08 mJy/arcmin 2 (right panel). For comparison, the size of the GBT field of 
view at this frequency is 32" x 32". 



in these pure forms, a realistic SZ flux decrement can be ob- 
tained by linearly combining the different scenarios due to the 
superposition property of the pressure of different populations: 

1 . The internal pressure is either dominated by CRps, mag- 
netic fields, or ultra-relativistic CRes being characterized 
by a mean momentum of (p) » 1 in this scenario. This 
is the most positive case for the detection of the plasma 
bubbles in the SZ flux decrement, as the bubble volume 
does not contribute to the SZ flux decrement significantly. 
Although this is the most positive scenario, it also is the one 
adopted by most analyses of radio lobes. 

2. This scenario assumes the internal pressure to be dominated 
by a compound of CRp and CRe populations where the lat- 
ter is described by a power-law distribution: 

fcR e (p,a,p u p 2 ) = —— 7—. (29) 

P\ ~P 2 

The distribution function is normalized such that its in- 
tegral over momentum space yields unity. The choice of 
Pi = 1 , P2 = 10 3 , and a = 2 implies a mean mo- 
mentum of (p) =i 6.9 as well as a pseudo-temperature of 
kfcRe — 2.2m e c 2 =i 1.1 MeV and represents a plausible 
scenario for the relativistic composition of the bubble. The 
CRp and CRe populations each contribute equally to the in- 
ternal pressure of the bubble representing a remnant plasma 
originating from the hadronic jet scenario. 

3. The dynamically dominant internal pressure support is 
contributed to equal amounts by relativistic electron and 



positron populations, respectively. Taking the same param- 
eters for the CRe distribution of the previous scenario, this 
approach represents the remnant radio plasma originating 
from the electron-positron jet scenario. 

4. A trans-relativistic thermal proton and electron distribution 
with kT e = 50 keV dominates dynamically over the other 
non-thermal components: 

/e,,h(p,Ah) = -^^P 1 eX P ("Ah V 1+ P 2 ) • (3°) 

K 2 denotes the modified Bessel function of the second 
kind (Abramowitz & Stegun 1965) which takes care of the 
proper normalization and f3 t h = m e c 2 /(kT e ) is the normal- 
ized thermal beta-parameter. The mean momentum of this 
distribution amounts to (p) 0.55. 

5. This scenario assumes a dynamically dominant hot thermal 
proton and electron distribution with kT s = 20 keV which 
exhibits a mean momentum of (p) =s 0.33. 

In all cases, the bubble's sound velocity and thus the baro- 
metric scale height are much higher than the corresponding val- 
ues of the ambient ICM which implies a flat pressure distribu- 
tion within the plasma bubble. This leads to a reduced internal 
pressure with respect to the ambient ICM at the inner rim of 
the bubble and an excess pressure at the outer rim of the bub- 
ble, which is responsible for the buoyant rise of the bubble in 
the cluster atmosphere. 

Figure 3 shows the spectral distortions due to the thermal 
SZ effect g(x), kinetic SZ effect h(x), and relativistic SZ ef- 
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Fig. 3. Spectral distortions due to the thermal SZ effect g(x), kinetic SZ 
effect h(x), relativistic SZ effect due to a population of ultra-relativistic 
CRes, gucRe = -i(x){3 VC R e (with p { = 3, p 2 = 10 3 , and a = 2), 
relativistic SZ effect due to a population of power-law CRes, g C Rc = 
[j( x )-i(x)]/3cRe (with/7[ = 1,P2 = 10 3 , and a = 2), and the relativistic 
SZ effect due to a population of trans-relativistic thermal electrons, 
§50 kcv = lj( x ) ~ /(x)]Ah(50keV), as well as due to electrons with 
kT e = 20 keV, respectively. 



feet of the various scenarios for the relativistic populations of 
the bubble composition. The spectral distortions of the rela- 
tivistic SZ effect are given by g (x) which have been defined in 
Eqns. (11) through (15). Please note, that the scenarios 2 and 3 
exhibit the same spectral distortion g cR e and that the amplitude 
of the kinetic SZ effect w gas is typically one order of magnitude 
smaller than the amplitude of the thermal SZ effect y gas . 

The left panel of Fig. 4 shows the unconvolved SZ flux 
decrement along an impact parameter through the center of the 
southern bubble of Perseus at the ALMA frequency band cen- 
tered on vo = 144 GHz. While the northern bubble of Perseus 
would evince qualitatively the same behavior, it shows a shal- 
lower depth of its SZ cavity due to its smaller geometrical size 
resulting in a weaker SZ flux contrast (cf. Fig. 1). The depth 
of the SZ cavity at this frequency range is a measure how rel- 
ativistic the respective electron population is, i.e. a deeper SZ 
cavity indicates a higher mean momentum of the electron pop- 
ulation. Our studies show a strong signature of the different 
bubble compositions on the SZ flux decrement. Thus, the com- 
bination of X-ray and SZ observations allows one to circum- 
vent the degeneracy between the effects of the bubble compo- 
sition and of the bubble extent along the line-of-sight on the SZ 
measurement. This enables us to distinguish a relativistic from 
a thermal electron population inside the bubble using only a 
single frequency SZ observation by either a detection or non- 
detection of the bubble, respectively. 



7. Kinetic Sunyaev-Zel'dovich effect 

7.1. General considerations 

If the cluster is moving towards the observer, the CMB tem- 
perature in this direction is increased due to the Doppler effect 
of the bulk motion of the cluster relative to the rest frame of 
the CMB (in this case, our convention is such that t> gas < 0). 
This enhanced SZ emission implies a reduced SZ flux decre- 
ment for frequencies below the crossover frequency of the 
thermal SZ effect at v c 217 GHz (for the non-relativistic 
case). Cosmologically, line-of-sight cluster velocities follow a 
Gaussian distribution with vanishing mean since the large scale 
structure is at rest in the comoving CMB-frame and with a stan- 
dard deviation of cr ve i ^ 310 km s~ as derived from a cosmo- 
logical structure formation simulation comprising the Hubble 
volume (Jenkins et al. 2001; Schafer et al. 2004a). 

Theoretically, the turbulent motions of the ICM should also 
contribute to the kinetic SZ effect. The largest impact is caused 
by two merging clusters along the line-of-sight which induces 
a bimodal streaming flow pattern on the sky. However, radio 
plasma bubbles are mainly observed within cool-core clusters 
which represent relaxed clusters where the ICM is approx- 
imately in hydrostatic equilibrium with the underlying dark 
matter potential. In this case, only small scale turbulent vor- 
tices are expected which exhibit smaller angular momenta. For 
an isotropic distribution of vortex orientations, line-of-sight in- 
tegration and beam convolution average the resulting kinetic 
SZ effect out to zero. Thus, these turbulent motions can only 
contribute in second order to the thermal SZ effect. 

It would be useful to know the peculiar velocity of the 
galaxy cluster in order to remove the degeneracy of the SZ 
cavity depth with respect to the different bubble compositions 
and the kinetic SZ effect. In the case of a general morphol- 
ogy of a galaxy cluster, in particular for non-axis-symmetric 
objects, it is impossible to unambiguously deproject the clus- 
ter in order to derive the peculiar velocity given only a single 
frequency SZ observation and an X-ray image (Zaroubi et al. 
1998, 2001). For multi-frequency SZ observations of the en- 
tire cluster, Aghanim et al. (2003) theoretically discuss possi- 
bilities in order to break parameter degeneracies between the 
Compton parameter, the electron temperature and the cluster 
peculiar velocity for an appropriate choice of observing fre- 
quencies. Considering possible gaseous substructure along the 
line-of-sight towards the radio plasma bubble and a possibly 
non-spherical bubble geometry, it might be impossible to dis- 
criminate between morphologically similar bubble composi- 
tions (scenarios 4 and 5) in the case of a single frequency ALMA 
or GBT observation. 

7.2. Perseus plasma bubbles at different frequencies 

The left panel of Fig. 4 compares three set of lines correspond- 
ing to different average line-of-sight streaming velocities at the 
ALMA frequency band which samples the extremum of the 
thermal SZ flux decrement (vo = 144 GHz). The effect of an 
overall decrease in the SZ flux decrement for the approaching 
cluster (y gas = -300 km s~') can be clearly seen. The ampli- 



C. Pfrommer, T. A. Enfilin, & C. L. Sarazin: Unveiling the composition of radio plasma bubbles 



9 



relative angular displacement [arcsec] 

20 40 60 80 



100 



ultra-relativistic CRe 

power-law CRe, e~- p - composition 

power-law CRe, e~- e + - composition . 

thermal electrons (kT c = 50 keV) 

thermal electrons (kT c = 20 keV) 




10 20 30 

relative distance [fc' kpc] 



40 



2.o r 



x 

o 



1.2 - 
1.2 - 



-2.0 L 



relative angular displacement [arcsec] 

20 40 60 80 



100 



= +300 km s" 




ultra-relativistic CRe 
power-law CRe, e~- p - composition 
power-law CRe, e~- e + - composition 
thermal electrons (kT c = 50 keV) 
thermal electrons (kT. = 20 keV) 




300 km s" 



10 20 30 

relative distance [/j 7( J kpc] 



40 



Fig. 4. Unconvolved SZ flux decrement along an impact parameter through the center of the southern bubble of Perseus. The left panel shows 
the SZ flux decrement at the ALMA frequency band centered on v = 144 GHz while the right panel assumes a central fiducial frequency of 
vo - 217 GHz. Compared are five different scenarios of the composition of the plasma bubbles to the undisturbed SZ profile (dotted thin line), 
respectively. The three (two) set of lines correspond to three (two) differently assumed average bulk velocities along the line-of-sight, y gas , of 
the thermal gas of Perseus. 



tude of the kinetic SZ effect depends on the dimensionless av- 
erage velocity of the thermal gas, which amounts in our case to 
ySgas = i> gas /c =s 1(T 3 . On the other hand, the amplitude of the 
thermal SZ effect depends on the inverse normalized thermal 
beta-parameter, f3^ = kT e /(m e c 2 ) 10~ 2 , where we inserted 
the average temperature of the Perseus cluster. At the extremum 
of the thermal SZ flux decrement, the kinetic SZ flux amounts 
to an approximately 10% correction to the thermal SZ effect for 
the choice of our fiducial cluster velocity. However, this small 
effect is responsible for the qualitative difference of the observ- 
able SZ cavities resulting from plasma bubbles: assuming a re- 
ceding cluster and a dynamically dominant thermal proton and 
electron distribution (scenarios 4 and 5), we still obtain a de- 
tectable depth of the SZ cavity which almost disappears for an 
approaching cluster. 

The right panel of Fig. 4 shows the unconvolved SZ flux 
decrement assuming a central frequency of vo = 217.34 GHz 
and bandwidth of Av = 40 GHz. Allowing for finite frequency 
response of the instrument's receivers, this fiducial frequency 
corresponds to a vanishing frequency band-averaged thermal 
SZ flux decrement. Thus, the kinetic SZ effect represents the 
main effect at this frequency showing a positive SZ flux decre- 



ment for a receding cluster, i.e. we would detect a reduced flux 
of CMB photons in that direction. At this frequency range, a 
relativistic electron population always causes a positive rela- 
tivistic SZ flux decrement owing to the higher crossover fre- 
quency v c > 217 GHz and irrespective of the cluster's velocity. 
This effect causes an enhanced SZ flux contrast of the bubble 
for an approaching cluster compared to a receding one, depend- 
ing on the specific composition of the plasma bubble. 

Most remarkable is the interchange of SZ flux decrements 
for different bubble compositions compared to the case of 
vo = 144 GHz. The largest impact of the bubble's SZ cavity is 
provided by the trans-relativistic thermal electron population of 
50 keV owing to its comparably large frequency band-averaged 
SZ flux decrement. The SZ flux decrements of the other bubble 
compositions are smaller because the crossover frequency of 
the thermal electron population of kT e = 20 keV lies closer to 
the non-relativistic crossover frequency while the amplitudes of 
the spectral distortions of the relativistic electron populations 
are smaller (cf. Fig. 3). Relativistic corrections to the thermal 
SZ flux resulting from the thermal electrons with temperatures 
of kT e ^3-7 keV are negligible: The resulting profiles of the 
SZ flux decrement are similar in morphology to the kinetic SZ 
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Fig. 5. Unconvolved SZ flux increment along an impact parameter 
through the center of the southern bubble of Perseus at the frequency 
band centered on vo — 470 GHz. Compared are five different scenar- 
ios of the composition of the plasma bubbles to the undisturbed SZ 
profile (dotted thin line), respectively. The two set of lines correspond 
to two differently assumed average bulk velocities along the line-of- 
sight, D g . ds , of the thermal gas of Perseus. 

flux decrement and would correspond at this frequency to an 
additional kinetic SZ effect of only ~ 10 km s . 

Owing to the proximity of the Perseus cluster to the 
Galactic plane, SZ observations at even higher frequencies 
might be challenging due to the Galactic dust emission which 
represents the major Galactic foreground at these frequencies 
(see e.g. Schlegel et al. 1998; Finkbeiner et al. 1999, 2000; 
Schafer et al. 2004b). In an exemplary manner, we also show 
the expected SZ cavity around a fiducial frequency which sam- 
ples the maximum of the SZ flux increment. 3 For a cluster 
with a temperature of 5 keV, the normal thermal relativistic 
corrections disappear at vo = 470 GHz when taking into ac- 
count the detector's finite frequency response of Av = 40 GHz. 
Secondly, at this frequency the kinetic SZ effect has only a 
small additional contribution (Aghanim et al. 2003). Hence, 
one could expect that the SZ signal from the bubble could 
be even more prominent. This is investigated in Fig. 5 which 
shows the unconvolved SZ flux increment along an impact pa- 
rameter through the center of the southern bubble of Perseus 
at vo = 470 GHz. The different bubble compositions qualita- 
tively show the same behavior as for the SZ decrement at the 
frequency band around vo = 144 GHz. However, the SZ flux 
contrast of the bubble is much higher as well as the degener- 
acy owing to the kinetic SZ effect is reduced. Thus, for clusters 
not being outshone by the Galactic dust emission, observations 

3 This might find application for a galaxy cluster with an X-ray cav- 
ity at high Galactic latitude (e.g., Abell 2052) or in the case of Perseus, 
once the Galactic dust emission is properly removed. 



at this frequency range seem to be most promising in order to 
infer the bubble composition. 

8. Observing strategy 

The previous results imply the following observing strategy for 
the ALMA compact core configuration or the GBT. A short du- 
ration observation (a few hours) of the inner rim of the south- 
ern X-ray cavity within the Perseus cool-core region should 
either yield a 5<x detection of an SZ cavity or not. The case 
of a significant detection either proves the existence of a dy- 
namically dominating CRp population, magnetic fields, or an 
ultra-relativistic CRe population within the radio plasma bub- 
ble while the contrary excludes this hypothesis on the 5cr level. 

Such a non-detection of the bubble would require a longer 
integration time in order to detect the SZ flux decrement of 
the plasma bubble at the desired significance level. A detection 
of a very shallow SZ cavity indicates a dynamically dominant 
hot thermal electron population within radio plasma bubbles. 
The underlying physical scenario is degenerate for an observed 
intermediate depth of the SZ cavity (such as an SZ flux level 
in between our scenarios 3 and 4). Possibilities include trans- 
relativistic thermal electron distributions or soft CRe power- 
law distributions exhibiting a spectral index of approximately 
a — 3 as well as a small lower cutoff p\ implying a mean mo- 
mentum of (p) ~ 1. Follow-up multi-frequency SZ observa- 
tions in combination with X-ray spectroscopy could disentan- 
gle the different scenarios and possibly estimate the temper- 
ature or spectral characteristics of the dynamically dominat- 
ing electron population (Shimon & Rephaeli 2002; Enfilin & 
Hansen 2004). The detection of such an SZ flux would enable 
one to draw conclusions concerning the particular jet scenario 
being responsible for the inflation of the plasma bubble. 

The case of a deep SZ cavity with an associated detec- 
tion of dynamically dominating CRps, magnetic fields, or ultra- 
relativistic CRes leads to an immediate question about the com- 
position of radio ghosts. Thus, an additional SZ observation of 
the ghost cavity in Perseus could yield answers about the po- 
tential entrainment of ICM into the plasma bubble during its 
buoyant rise in the cluster atmosphere. If a large fraction of 
entrained gas in the bubble provides significant pressure sup- 
port, it must have experienced Coulomb heating through CRes 
in order not to be detected in X-rays. This would result in a re- 
duced SZ flux contrast of the ghost cavity and leads to a faint 
or undetectable SZ flux decrement of the ghost cavity. 

9. Conclusion and outlook 

This work provides a theoretical framework for studying the SZ 
decrement of radio plasma bubbles within clusters of galaxies. 
X-ray observations are proportional to the square of the ther- 
mal electron density and probe the core region of the cluster. 
On the other hand, the SZ effect is proportional to the thermal 
electron pressure enabling the detection of plasma bubbles fur- 
ther outwards of the respective clusters. Assuming spherically 
symmetric plasma bubbles, we simulate an ALMA and a GBT 
observation of the cool core regions of the Perseus cluster and 
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Abell 2052. In this context, we exemplarily investigate physi- 
cally different scenarios of the composition of the plasma bub- 
bles: As long as the bubble is dynamically dominated by rela- 
tivistic electrons, protons, or magnetic fields, there exists a real- 
istic chance to detect plasma bubbles in the SZ flux decrement. 
Non-detection of radio bubbles with the SZ effect at the po- 
sition of X-ray cavities hints towards a dynamically dominant 
hot thermal electron population within radio plasma bubbles. 
Detection of a non-thermal pressure support should be possible 
within a few hours observation with ALMA or GBT in the case 
of Perseus and a few ten hours in the case of Abell 2052. 

Pursuing high-sensitivity multi-frequency SZ observations, 
it will be challenging but not impossible to infer the detailed 
nature of the different possible populations of the bubble. For 
realistic observations, the frequency dependence of the rela- 
tivistic SZ signal is contaminated by the presence of the kinetic 
SZ signal. It would be optimal to have a frequency channel 
centered on the crossover frequency of the thermal SZ effect 
in order to infer the energy spectrum of the dynamically dom- 
inating population and to measure the temperature or spectral 
characteristics of the electron population. 
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Appendix A: Synthetic ALMA observation of the 
central region of Perseus 

This appendix provides additional details of the synthetic 
ALMA observation of the cool core region of Perseus as well 
as the expected SZ flux contrast of the southern X-ray cavity. 
First, we compute the frequency band-averaged SZ flux decre- 
ment which is defined by 



/ dvSI(v)R Vo (v) 
fdvR Vo (v) 



(A.l) 



Here, R Vo (v) denotes the frequency response of the ALMA re- 
ceivers centered on a fiducial frequency vo which we assume to 
be described by a top-hat function: 



#V (V) = 



l,ve [v - Av/2, v + Av/2] 
0, v $ [v - Av/2, v + Av/2] 



(A.2) 



We choose the ALMA frequency band 4 which samples the 
extremum of the SZ flux decrement and is characterized by 
vo = 144 GHz and Av = 38 GHz (Brown et al. 2000). The 
requirement of obtaining the highest flux sensitivity to the 
largest scales comparable to the field of view at this frequency 
(FWHMf ov =2 36") calls for the most compact configuration 
ALMA E with a maximal baseline of <4i = 150 m. Thus, we 
convolve the simulated SZ flux decrement with a Gaussian 
to obtain the resolution of this configuration, FWHM 
c/(v dbi) - 2.86". Assuming an ultra-relativistic electron popu- 
lation within the bubbles yields the SZ flux decrement as shown 
in Fig. 1. 



To investigate, whether the plasma bubbles are detectable 
in the SZ flux decrement, we define the SZ flux contrast p: 



P = 



Ia-Ib 



(A3) 



Here, I A and / B denote the mean SZ flux decrement of two 
equally sized solid angle elements within the field of view, one 
of which measures the SZ flux inside and the other one out- 
side the bubble. Adopting Gaussian error propagation and in- 
troducing ALMA's sensitivity per beam in terms of flux density, 
(Talma, yields the uncertainty in p, 



(A.4) 




2 fj 2 + J 2 
\ a b ctalma 



(7 a + /b) 2 VA^, 



(A.5) 



Here, A^eam = /" Abeam /(2Af ov ) is the number of statistically 
independent beams per flux averaged solid angle element, / < 
0.5 measures the fraction of the solid angle of the bubble within 
the field of view, and the factor of Abeam /(2Af ov ) accounts for 
the statistically independent degrees of freedom within the field 
of view. Combining Eqns. A. 3 and A.5 yields the signal-to- 
noise for the 5cr detection of the bubble: 
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(A.6) 



The flux density sensitivity (Talma for point sources, which 
should be approximately applicable in the case of the compact 
core configuration, is given by Butler et al. (1999): 



calma - 



V2>r sys 



^Adish A beam VAf AvA^bi 



(A.7) 



Here, r sys ^ 65 K denotes the system temperature, Ad; s h = 
7rD 2 /4 the collecting area of each dish (D = 12 m), rj = 0.75 
the aperture efficiency, Abeam = 1.13 FWHM 2 9.3 arcsec 2 
the area of the secondary beam, Af the integration time, Av = 
38 GHz the bandwidth, and A^i = n(n - l)/2 the number of 
baselines, n = 64 being the number of antennas. We choose 
the field of view to be centered on the inner rim of the south- 
ern bubble such that equal solid angle elements fall inside and 
outside the bubble, respectively. The simulated mean SZ flux 
decrements within the two solid angle elements of the primary 
beam, I A = 13.25 mJy arcmkT 2 and 7b = 12.70 mJy arcmin -2 , 
require only an integration time of 5.1 hours in order to ob- 
tain a 5cr detection of the plasma bubble assuming an ultra- 
relativistic electron population within the bubble. Uncertainties 
in the amplitude of the kinetic SZ effect and the geometrical ar- 
rangement and shape of the bubbles may slightly modify this 
result. 
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